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Abstract 

We analyze two-loop gluino corrections to the neutron electric dipole moment (EDM) in the 
minimal supersymmetry extension of the standard model (MSSM). The dependence of two- loop 
corrections on the relevant CP violating phases differs from that of the one- loop contributions, and 
there is a region in the parameter space where the two-loop contributions are comparable with the 
one-loop contributions. Our numerical results show that the two-loop corrections can be as large 
as 30% of the one- loop results. 
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I. INTRODUCTION 



The fermion electric dipole moments (EDMs) offer a powerful probe for new physics 
beyond the Standard Model (SM). In the SM, the EDM of the neutron is fully induced by the 
CP phase of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and it is predicted 
to be much smaller [l] than the present experimental upper limit of 1.1 x 1CT 25 e ■ cm Q] 
and beyond the reach of experiment in the near future. As for the minimal supersymmetric 
extension of the SM (MSSM), there are many new sources of the CP violation that can 

nn 

result in larger contributions to the EDM of the neutron jJJ ]4^. Taking the CP phases 

with a natural size of 0(1), and the supersymmetry mass spectra at the TeV range, the 

theoretical prediction on the neutron EDM at one-loop level already exceeds the present 

experimental upper bound. In order to make the theoretical prediction consistent with the 

experimental data, three approaches are adopted in the literature. One possibility is to make 

pi 

the CP phases sufficiently small, i.e. < 10~ 2 |2j]. One can also assume a mass suppression 
by making the supersymmetry spectra heavy, i.e. in the several TeV range or invoke a 
cancellation among the different contributions to the fermion EDMs j^J . 

The prediction for the fermion EDMs at one-loop level in a supersymmetric theory has 
been extensively discussed in the literature. On other hand, the analysis on the neutron 
EDM at two-loop level is less advanced, even though some pioneer work has been carried 
out, for example the two-loop Barr-Zee type diagrams involving the Higgs bosons p] and the 
purely gluonic dimension-six Weinberg operator induced by two-loop gluino-squark diagrams 
j7( have been analyzed. 

Analyzing the EDMs at two-loop order can give us a better understanding of where the 
new physics scale may emerge and shed some light to the spectra of new physical particles 
around this scale. Moreover, the two-loop analysis involves some new parameters in addition 
to those appearing in the one-loop calculations, hence resulting in more rigorous constraints 
on the supersymmetry parameter space. 

In this work we shall analyse the two-loop gluino corrections to the neutron EDM in the 
MSSM. We work in the framework of the simplest model, where we will neglect all other 
possible sources of flavor violation except those related to the CKM matrix, and try to 
avoid ambiguities of the unification conditions of the soft-breaking parameters at the grand 
unification scale present in the mSUGRA schemes. In the next Section, we will demonstrate 
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how to obtain the two-loop gluino corrections to the neutron EDM. In the Section IIHI we 
will present the results of our numerical computations and we study the dependence of the 
neutron EDM on the supersymmetry parameters. Conclusions are presented in Section IIV1 

II. THE TWO-LOOP GLUINO CORRECTIONS TO THE NEUTRON EDM 

In the effective Lagrangian, the fermion EDM d, is defined through the dimension five 
operator 

£ EBM = ~d f Ja^ l5 fF^ (1) 

where/ is a fermion field, F is the electromagnetic field strength. This coupling obviously 
signifies a CP violation. It is not present among the fundamental interactions at tree level, 
but it is generated in loop level in an electroweak theory with CP violation. Moreover, 
because quarks also take part in strong interactions, the chromoelectric dipole moment 
(CEDM) fT a a^ u ^^fG a ^ of quarks, where T a (a = 1, • • • , 8) denote the generators of the 
strong SU(3) gauge group, G a v is the gluon field strength, and the pure gluon Weinberg 
operator of dimension- six, f „ G a G bp G c e^ A<J , contributes to the quark EDMs as well. 

A convenient way to describe loop-induced contributions is the effective theory approach, 
where the heavy particles are integrated out at the matching scale. The resulting effective 
Lagrangian includes a full set of CP violation operators. In this work, we restrict ourselves 
to the following operators that are relevant to the neutron EDM: 

^, = £0^)0, (A), (2) 

i 

where C.(A) are the Wilson coefficients evaluated at the scale A, and the five operators of 
interests are 

Qi = qa^qF^ , 
Q 2 = qa^co+qF^ , 
Q 3 = qTV^qGl, 
Q 4 = qWu; + qGl , 

2 5 = ~h b GlG b fGle^ . (3) 
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FIG. 1: The one-loop self energy diagrams which lead to the quark EDMs and CEDMs in MSSM, 
the corresponding triangle diagrams are obtained by attaching a photon or a gluon in all possible 
ways to the SUSY particles. 



The effective triangle diagrams responsible for the quark EDMs and CEDMs, are obtained 
by attaching a photon or gluon external line,for quark EDM or CEDM, respectively, in all 
possible ways to the quark self-energy diagrams. 

The one-loop supersymmetric corrections to the Wilson coefficients Cj(A) in Eq. (J2J) 
originate from three types of the graphs: the quark self-energy diagrams, where one has 
a gluino and a squark, a chargino and a squark or a neutralino and a squark as virtual 
particles in the loop (FigOJ. We will adopt below a terminology where, for example, the 
"gluino- squark contribution" means the sum of those triangle diagrams, which have a gluino 
and a squark as virtual states and where a photon or a gluon in attached in all possible ways 
to the external line. 

Before moving to derive the two- loop gluino corrections, we first present, in order to fa- 
miliarize our notation, the one-loop results. Those expressions can be found in the literature, 
and we just translate them here into our notations for later convenience. 

The corrections from the gluino-squark diagrams are 

2 2 

d? = e ea V Jm((ZX XZt).,e 

° 11 1=1 
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m 2 v m 



Here a g = g 2 /(ATr), 6 :j denotes the phase of the soft gluino mass m., and Z. are the mixing 
matrices of the squarks, i.e. Z\m 2 Z~ = diag(m 2 , m 2 ) where 
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m 2 

9 



m 2 +m 2 + m 2 (i-Q s 2 ) cos 2/3 m (A* - uR ) \ 

m q (A q -n*R q ) va\. b) + w? + m 2 (± - Q qS 2 ) cos 2/3 ) 



, (5) 



{£/,£>} 

with Q q = 2/3(— 1/3), i? = tan/3(l/ tan/3) for g = w (d). As usual, tan/3 = v u /v d is 
the ratio between the VEVs of the up- and down-type Higgs fields, and 9 W is the weak 
mixing angle. We also use the short-hand notations s w = sin6> w , c w = cos# w . The loop 
functions B(r) and C(r) are defined as B(r) = [2(r — 1) 2 ] _1 [1 + r + 2rlnr/(r — 1)] and 
C{r) = [6(r - l) 2 ]" 1 ^^ - 26 - (2r - 18) lnr/(r - 1)]. 

In a similar way, the one-loop neutralino-squark contributions can be written as 

m x k , m2 k 



with 



2 \ m 2 / ' 



x- , 



d9 = 93a ° Ylm((A q ) h .(B q y 
*° k m 64vrs 2 c 2 4? " NJkA Nh '< 



m , m 2 



^BO («) 



4 772 C 

(^")fc,i = _ o S w(^'ii)2.i(^'jv)l.fc + ZT" 

m w 6 /3 

X (^)i,i(^Ar) 4 ,fe ) 
(^)w = (^i.«(T( Z X + C -( Z X) 



(^)fc.i _ q S w(^ d -) 2 , l (^ : iv)i,fc + " 

m w C /3 

X(^)i. ( (^)s lfc , 

(^)M = (^) 1 , l (f(^):, fc -c w (^):j 
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Here a = e 2 /(47r), s = sin/3, c„ = cos/3, m k (k — 1, 2, 3, 4) denote the eigenvalues of 

x o 

neutralino mass matrix, and Z N is the correspondingly mixing matrix. 
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Finally, the chargino-squark contributions are given by 
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Airs 2 
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Qi 



m 



m 



x 



Qi Qi 
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q Q Qq C'k,i\ c'i,k 



93a 

Ans 2 

w 

m , 



x ■ 



m 



(8) 

Qi Qi 

where V denotes the CKM matrix, m , (k — 1, 2) are the masses of charginos, and the 
loop function A(r) = 2(1 — r)~ 2 [3 — r + 2 lnr/(l — r)]. The couplings appearing in these 
expressions are defined as 



y/2m. .s 



w"/3 



V2m w c p 



V^ d -)l,i(^+) 2 ,fe > 

(^) w (^-) a . fc -(^)l. 4 (^-)l>, 




(f) 

FIG. 2: The two- loop self energy diagrams which lead to the quark EDMs and CEDMs in the 
MSSM, the corresponding triangle diagrams are obtained by attaching a photon or gluon in all 
possible ways to the SUSY particles. 
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{z,)uz + ): 



(9) 



where Z+ and Z_ are the right- and left-handed mixing matrices of the charginos, respec- 
tively. It is noteworthy that one-loop chargino corrections to the quark EDMs and CEDMs 
are proportional to the suppression factors m q /m vi contained in the couplings A q c (q = u, d). 

The two-loop gluino corrections to the Wilson coefficients originate from the two-loop self- 
energy diagrams for quarks, depicted in Fig. 121 The corresponding dipole moment diagrams 
are obtained by attaching a photon or gluon to these diagrams in all possible ways. In these 
diagrams there is no new suppression factor, except a factor arising from loop integration, 



and the divergence caused by the sub-diagrams is subtracted in MS scheme safely. It turns 
out that for some region of the parameter space the two-loop results are comparable with the 
one-loop contributions. The reason for this is that the dependence of the two-loop results 
on the relevant CP phases differs from that of the one-loop results. 

Since the two-loop analysis is more subtle than the analysis at one-loop level, we present 
here in some detail all the processes, which contribute at two-loop level into the theoretical 
)rediction of the quark EDMs. Taking the same steps, which we did in our earlier works 
we obtain the following expression for the relevant effective Lagrangian: 
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+(^) 2 ,(^), 2 (2 f ) 2 , i (^) i , 1 |m § |e^3 
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e q e^ {2) Q 2+ g 3 W (2) Q 4 }} + 



(10) 



where T>. 



((92 - qi) 2 ~ m 2 )(qj - |mJ 2 )(g? - m 2 )(g| - m 2 _)(g 2 2 



\m. 

i a 



I 2 ) and D is the 



time-space dimension. We collect the complicated form factor M- in appendix EI 

The effective Lagrangian of Eq. (jTUj) originates in fact from the contributions of the 
effective triangle vertices induced by the "gluino-gluino self-energy" (Fig. 2 (a)). The termi- 
nology here is analogous to the one adopted in the one-loop case above. The gluino-gluino 
self-energy is a short-hand expression for a diagram, where we attach a photon or gluon line 



in all possible ways to the two-loop quark self-energy diagram having, besides the squark, 
two gluinos are as intermediate agents. Because the sum of the triangle diagrams corre- 
sponding to each "self-energy" obviously respects the gauge invariance, we can calculate the 
contributions of all the "self-energies" separately. We will use below the following identities 
demanded by the translational invariance of the loop-momenta: 

d D q 1 d D q 2 1 r 2 + D 
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where V = ((q 2 - qi) 2 - m 2 ){q\ - m 2 )(q 2 - m 2 ). We find 
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4 qfqi • (?2 - gi) , qi • (ga - 9i) , 9i • (?2 - ?i) 
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Similarly, we find 
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d°gl d% rnl 

(2tt) d (2tt)^ P s " 1 J 

With the help of Eq.(fT2"j) and Eq. tfTSj) . the corresponding terms in Eq.fITU]) are transformed 
into 

iy- 9 = _ % i e 4 |m | f dPqi d% ^ , .(^t) e -* 3 ) . [^75^' 

4 /• /7-°/r rl D n Af 9 

-i-gl\m„\ / ° gl n , f n ^ Imf(2:-) 2 (Zt).^-^) . \qT a a^ 5 qG a uu 
9 3 9 i (2tt) d (2tt) d V ^ q J 9 J ' ' L m J 

+ ■••• 9 (14) 

The two-loop scalar integrations can be reduced into the two-loop vacuum integrals defined 

as 

7_C> jD 



d D qi d u q 2 



(qi ~ g 2 ) 2 - rnl){ql - mf)(g| - m 2 ,) ' 

3D D AT 7 ' 9 

. However, the two- loop scalar integration / —j^ 1 

contain the ultra-violet divergence which originates from the sub-daigrams. After the renor- 



malization in MS scheme, we finally obtain the corrections from the two-loop " gluino-gluino" 
diagram to the quark EDMs and CEDMs respectively as 



(56 &0?\Tfl 1 
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J 3 
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with Xj = m 2 /m 2 , and the function -Fj(x , £1, x 2 , £3, x 4 ) (i = 1, 3) are defined in appendix 

m 

The effective Lagrangian, which corresponds to the contributions of the triangle dia- 
grams induced in the "neutralino-gluino" self-energy (i.e. diagrams where loop squarks are 
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accompanied with neutralinos and gluinos, see Fig.2 (b) and (c) ), includes the following 
pieces: 
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X°(l,2) 

With a help of Eq. (fTTj). one finds the following identities: 



d D qi 


d D q 2 


x°(i) x °m 


(2tt) d 


(2ir) D 


x° 
x k 


d D Ql 


d D q 2 


A/" 7 ,Af 9(a) 

9(2)' X°(2) 


(2tt) d 


(2tt) d 


£>(«) 

4 


d D gi 


d D q 2 


x°(i)' x°(i) 


(2n) D 


(2ir) D 


VQ>) 

< 

k 


d D qi 


d D q 2 


A/" 7 ,A/" 9(a) 

9( 2 )' X°(2) 


(2tt) d 


(2tt) d 


£>( a ) 


d D qi 


d D q 2 


A/" 7 ^ ,JV* 6 > 

X°(2)' x°(2) 



(ft 



ft, ft 



ft 



(ft 



gi <-> g2 



(2tt) d (2tt) 



D 



fe 



ft 



ft, ft «-> ft 



10 



d D Ql 


d D q 2 


A/ 7 ,A/" s(a) 

§(!)' X°(l) 


(2tt) d 


(2tt) d 




d D gi 


d D q 2 






(2tt)» 




d D Ql 


d D q 2 


A/" 7 M 9{a) 

§(!)' X°(l) 


(2ir) D 


(2tt) d 


*2 



(17) 



Substituting these identities into eq.(15) and removing the ultra-violet divergence in MS 
scheme, we may extract the two- loop "neutralino-gluino" corrections to the quark EDMs 
and CEDMs: 
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In a similar way, we get the following results for the two-loop "gluino-chargino" correc- 
tions: 

4eaa, 



d J = 
x±(2) 3(47r) 2 s 2 m 



d 9 , 



^^{l^l J? 6(x Q ,^, a ; B ,x^, a ; fl( )[lm((^) w (2j) Jil 

-m^F 7 (x g ,x aj ,x B ,^,x ft )[lm((A«) w (^) Jia (^) H (Zt) 4i ) 
-Im((S?) w (2j) J , 1 (S«) H (Zt) ia )]} > 



x*(2) 3(47r) 2 s 2 m 

x^'U^e"*") _ Im((S«) w (^) Jia (^) H (^) iiae «3)] 

^/5(^.^/,.^.^)[ Im ((4( 2 S) A .(^( 2 !)M) 



11 



-Im((^) fcj (Zt) j , 1 (^) fei (Zt) M )]}. 

(19) 



Notice that the last terms of the d 1 and d 9 are not proportional to the suppression 
factor m q /m w . This implies that the two-loop "gluino-chargino" diagrams may be dominant 
two-loop corrections of the chargino to the quark EDMs and CEDMs. 

The Wilson coefficient of the purely gluonic Weinberg operator originates from the two- 
loop "gluino-squark" diagrams and is given by jiol ]: 



G 6 = - 3 W£)M^U^)L 



\4 n J 
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\m. r \m.r \m~ r \m~ r 

I g I Iffl I9I Ifll 

where Z- t is the diagonalizing matrix for the squared mass matrix of stop, and the function 
H is the same as that given in [loj. 

In order to account for the resummation of logarithmic corrections, we should evolve 
the quark EDMs, CEDMs and the Wilson coefficients of the Weinberg operator with the 
renormalization group equations (RGEs) from the matching scale A down to the chirality 
breaking scale A x |ll| : 

C 5 (A x )= Vg C 5 (A), (21) 

where r] i ~ 1.53 and r] g ~ r] G ~ 3.4. At a low scale, the quark EDM consists of a weak 
interaction contribution, a quark CEDMs contribution and the contribution of the Weinberg 
operator. We can numerically evaluate these different contributions. According to a naive 
dimensional analysis one can write [12| 

e <°A 

d =d? + — d 9 + — ^CAA ) . (22) 

The EDM of the neutron is, on the basis of the simple SU(6) quark model, then given by 

d n = ±{4d d -dj, (23) 

where d u , d d are the EDMs of the up- and down-type quarks respectively. 
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FIG. 3: The correction from one-loop gluino diagrams (Eq. and that from two-loop gluino- 
gluino diagrams (Eq. 115(1 versus the right handed squark masses = with 9 3 = arg(m.) = 
7r/2, tan (3 = 6 (60), and m- = 20 TeV, where the dash-lines represent one-loop results and the 
solid-lines correspond to a sum of the one- and two-loop results. 

III. THE NUMERICAL ANALYSIS 

The MSSM Lagrangian contains several sources for CP violating phases: the phases 
of the /i parameter in the superpotential and the corresponding bilinear coupling of the soft 
breaking terms, three phases of the gaugino mass terms, and the phases of the trilinear 
sfermion Yukawa couplings in the soft Lagrangian. As we are not considering the sponta- 
neous CP violation in this work, the CP phase of soft bilinear coupling vanishes due to the 
tree level neutral Higgs tadpole conditions. 

As we remarked above, the two-loop gluino corrections to the neutron EDM are im- 
portant, and may be the dominant contribution in a certain part of the parameter space. 
Let us clarify this point in some detail by means of numerical analysis. Without losing 
too much generality, we will fix the following values for the supersymmetric parameters: 
M = l m A B l = l m A A l = \ m - g \ = 300 GeV, 1 \A q \ = 100 GeV (q = u, d, c, s, t, b), and 

1 m A , to a are the soft masses of U(l) x SU(2) gauginos respectively. 
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m. = 20 TeV, m~ =5 TeV, m~ — m~ =5 TeV, m. = m~ = 200 GeV in our numerical 

Q 2 ' Q 3 'OS ' T B 

analysis. For simplification, we will also assume (f q = arg(A q ) = 0, 8^ = arg(fi) = 0. 

In Fig. El we present the corrections the neutron EDM obtains from the one- loop "gluino- 
squark" diagrams (see Eq. (0J) and Fig. 1) and two-loop "gluino-gluino" contributions (see 
Eq. (JT5j) and Fig. 2) as a function of the right-handed squark masses m~ = . We have set 
in this figure 8 1 = arg(m XB ) =0, 9 2 — arg(m XA ) =0, 9 3 — arg{m g ) = 7r/2, tan/3 = 6 (60), 
and = 20 TeV, and what are plotted are the one-loop "gluino- squark" contribution 
(dashed line) and the sum of the one- loop " gluino-squark" and two- loop "gluino-gluino" 
contributions (solid line). The dependence of the contributions on the CP violation phase 8 3 
are all proportional to the factor Im^Z.) 2 ^(Z^) - 1 e~ l6 '^ , and the two-loop correction is less 
than 10% of the one-loop results for our choice of the parameter values. The neutron EDM 
depends on the parameter tan /3 through the squark mixing matrices, and as seen from the 
plots, this makes the a large difference between the corrections in a low and a high tan/3 
cases. 

Let us now focus on the two-loop gluino corrections of the neutron EDM. There the CP 
violating phases 8 1 , 8 2 induce a nonzero quark EDMs at one-loop level (see Eq. ©) through 
the neutralino mixing matrix, whereas the CP phase 8 3 contributes only through the two- 
loop diagrams (see Eq. (fTHjl ). We will present three sets of plots, where we set one of the 
phases 81,82,83 at a time equal to 7r/2 while the other two are set equal to zero. Taking 
first 6 1 = 7r/2, 8 2 = 8 3 = and tan/3 = 6 (60), we plot in Fig. H]for two values of tan/3 the 
neutron EDM as a function of the soft squark masses m. = m. = m.. In this plot the 

^ Ql U D r 

dash-lines represent the one-loop " neutralino-squark" results and the solid-lines correspond 
to a sum of one-loop "neutralino-squark" and two-loop " neutralino-gluino" contributions. 

In Fig. Elwe present the same plot for the case 8 2 = tt/2, 8 1 = 8 3 = 0, other parameters 
being the same as in Fig. |3J We notice that the two-loop gluino corrections to the one-loop 
results can be as large as about 20%. For the contribution of neutralino to the neutron 
EDM, the most interesting piece originates from the CP violating phase of the gluino mass, 
which does not contribute at the one-loop level. 

In Fig. El we plot the neutron EDM versus the soft squark masses = = m b in the 
case 8 3 = n/2, 9 1 = 8 2 = 0, the other parameters are taken as in Fig. |3J As the soft squark 
masses increase, the neutron EDM d n varies steeply. Because of the interference between 
the up- and down-quark EDMs (see Eq. I2l?|) . there is in the case of tan/3 = 6 a profound 



14 



1 E-28 I I 

1 2 3 4 5 

m g — m f. = m p (TeV) 

FIG. 4: The correction from one-loop neutralino diagrams (Eq. EJ) and that from two-loop 
neutralino-gluino diagrams (Eq. I18j) versus the squark masses = = with 9 ± = 

7r/2, 9 2 = 8 3 = 0, tan/3 = 6 (60), where the dash-lines represent one-loop results and the solid- 
lines correspond to a sum of the one- and two-loop results. 

minimum at m, = m . = m ~ ~ 900 GeV . 

Ql U D 

Let us move to the corrections arising from the diagrams involving virtual charginos. In 
this case the phase 6 X does not contribute to the quark EDMs at all, whereas 9 2 contributes 
at one- loop level (see Eq. (jSJ)) and 9 3 takes part in the game through only two- loop diagrams 
(see Eq. ()19)) ). At the one- loop level, the chargino contributions to the quark EDMs are 
proportional to the suppression factor m q /m z (q = u, d). However, the last terms of the 
two- loop "gluino-chargino" corrections to the quark EDMs and CEDMs, given in Eq. (|T^|) . 
do not depend on this suppression factor. Taking 9 2 = n/2, 9 3 = 0, we plot in Fig. 
the neutron EDM as a function of the soft squark masses m . — m~ = m, for a large 

H Ql V D to 

tan/3. The dash lines represent the one-loop "chargino- squark" results and the solid lines 
correspond to the sum of the one- loop "chargino- squark" and two-loop " chargino-gluino" 
results. As the dominant two- loop "gluino-chargino" corrections are not proportional to the 
suppression factor given above, they play the leading role. The corresponding figure for a 
low tan (3 is presented in Fig. |H1 
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1 E-28 I I 

1 2 3 4 5 

m g — m f. = m p (TeV) 

FIG. 5: The correction from one-loop neutralino diagrams (Eq. IHJ only and the sum of one-loop 
neutralino and two-loop neutralino-gluino diagrams (Eq. I18|) versus the squark masses m ~ = 
m & = with 9 2 = vr/2, 9 1 = 3 = 0, tan/? = 6 (60), where the dash-lines represent one-loop 
results and the solid-lines correspond to a sum of the one- and two-loop results. 



Just as in the neutralino case, the CP phase 9 3 affects also in the chargino case the quark 
EDMs only at the two-loop level. In Fig. El we show the variation of the neutron EDM with 
respect to = m = m 6 in the case 9 3 = it/ 2, 9 2 = and tan/5 = 6 (60). Even when 
the squark masses are larger than 1 TeV, the neutron EDM induced by the CP phase 9 3 can 
be close to 1.1 x 10 -25 (e • cm), which is the present upper bound from experiments. 

By including all the contributions present in Eq. (|22jh we now determine in our theoretical 
framework the variation of the neutron EDM as a function of the CP phase 9 3 . Taking 
6 1 = 9 2 = 0, = = m fi = 2.3 TeV and tan/3 = 60, we plot in Fig. El the neutron 
EDM versus the CP phase 9 3 , where the dashed line is the sum of the one-loop results 
and the contributions of the Weinberg operator, and the solid line represents the results 
including two-loop gluino corrections to the quark dipole moment operators. Fig. ITT1 gives 
the corresponding plot for the case tan (3 = 6. It is clear that the two-loop gluino corrections 
can reach the 30% level for our choice of the parameters. 

In this work, we do not consider the effects of the CP phases associated with the trilinear 
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12 3 4 

m g — m f. = m f, (TeV) 

FIG. 6: The correction the two-loop neutralino-gluino diagrams (Eq. I18|) versus the squark masses 
= = with 9 3 = vr/2, 9 1 = 9 2 = 0, tan/3 = 6 (60). With the choice on parameter 
space, the one-loop neutralino contribution to neutron EDM is zero. 




1 E-27 I I 

1 2 3 4 5 

TTi - = rn,-_ = (TeV) 

FIG. 7: The correction from one- loop chargino diagrams (Eq. |HJ) and the sum of one- and two- loop 
chargino-gluino diagrams (Eq. I19|) versus the squark masses = = with 6 2 = ir/2, 8 3 = 
0, tan/? = 60. 
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1 E-27 I I 

1 2 3 4 5 

m g = m. = (TeV) 

FIG. 8: The correction from one-loop chargino diagrams (Eq. EJ) and the sum of one- and two-loop 
chargino-gluino diagrams (Eq. PHI)) versus the squark masses = = with 9 2 = tt/2, 8 3 = 
0, tan/? = 6. 




FIG. 9: The correction from the two-loop chargino-gluino diagrams (Eq. I19|) versus the squark 
masses m- = m~ = m ~ with 9 3 = ir/2, 1 = 2 = 0, tan/? = 6, (60). 
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0.0 0.5 1.0 1.5 2.0 

FIG. 10: Taking 1 = 6 2 = 0, = m. = m. = 2.3 TeV and tan /? = 60, the neutron EDM 

varies with the CP phase 6 3 . 
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-10 i i 

0.0 0.5 1.0 1.5 2.0 



FIG. 11: Taking 1 =0 2 = 0, m. = m. = m. = 2.3 TeV and tan/3 = 6, the neutron EDM varies 
with the CP phase 6 3 . 
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soft squark Yukawa coupling. This is reflected in our numerical results as a relatively strong 
decrease of the neutron EDM as a function increasing squark masses. Moreover, for the 
model we employ here, the mass of the lightest Higgs boson sets a strong constraint on the 
parameter space of the new physics. As indicated in the literature [13j, the CP violation 
would cause changes to the neutral-Higgs-quark coupling, neutral Higgs-gauge-boson cou- 
pling and self-coupling of Higgs boson. The present experimental lower bound for the mass 
of the lightest Higgs boson is relaxed to 60 GeV. In our numerical analysis we have taken 
this constraint for the parameter space into account. 

IV. THE CONCLUSION 

We have investigated in this work the corrections induced by some two-loop gluino 
diagrams to the neutron EDM. Except a loop factor, the two-loop contributions do not 
involve other suppression factors. Since the dependence of the two-loop corrections on the 
relevant CP violating phases differs from that of the one-loop corrections, there is a region 
in the parameter space, where the two-loop corrections are comparable with the one-loop 
corrections. Certainly, the present experimental upper bound on the neutron EDM would set 
a very rigorous constraint on the possible values of CP phases if the masses of new particles 
are in the TeV range. In order to circumvent this constraint, one can assume, as was done 
in the corresponding analysis at one-loop level, that the squarks of the first two generations 
are heavy, i.e. > 2 TeV, or invoke a mechanism that causes an effective cancellation among 
different contributions to the quark EDMs. 

Because the exact analytical expressions for the two-loop vacuum integrals are derivable, 
one can extract, at least in principle, the Wilson coefficients of the dipole moment operators 
from any two-loop triangle diagrams and thereby find their contributions to the fermion 
anomalous magnetic dipole moments and EDMs. We will address our analysis to this issue 
in our coming works. 
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APPENDIX A: THE FORM FACTORS 
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APPENDIX B: THE EXPRESSION OF F i (x ,x 1 ,x 2 ,x 3 ,x 4 ) (i = 1, 2, • • • , 7) 



Fi(xo,xi,x 2 ,X3,x 4 ) = {S A + E B + E C + E D + E E + is F }(x , xi, x 2 , x 3 , x 4 ) , 

3 1 

F 2 (x ,Xi,x 2 ,a;3,X4) = {E A + -S s + E c + S D + S B + -H F }(x , x 2 , xi, x 4 , x 3 ) 

F 3 (x ,x 1 ,x 2 ,x 3 ,x 4 ) = {E a + yS s + S C - + E E - ^E F }(x ,x 1 ,x 2 ,x 3 ,x 4 ) 

9 r 

-Jq\£ b (x ,X 2 ,X 1 ,X4,X 3 ) + E B (x ,X4,Xi,X 2 ,X 3 ) 

F 4 (xo,Xi,x 2 ,X3,x 4 ) = {S A + S B + S c - -S c + E E - —E F }(x ,x 1 ,x 2 ,x 3 ,x A ) , 
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F 5 (x ,xi,x 2 ,x 3 ,x 4 ) = {E A + + S c + E D + E E - ^S F }(x , x 2 , x u x 4 , x 3 ) 
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Where the functions S a (xo, x\, x 2 , x 3 , x 4 ) (a = A, B, C, D, E, F) are formulated as 
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Z F (x ,x 1 ,x 2 ,x 3 ,x 4 ) = ^{ ^ _ - ^ _ - j $(x ,x 1 ,x 2 ) - $(x ,Xi,X 4 ) 

Xo + x\ — x 2 (9$ 



-$(x , x 3 , x 2 ) + $(x , x 3 , x 4 ) 
xq + xi — ^4 9$ 



+ 



(xi - x 3 )(x 2 - x 4 ) dx 



(x ,xi,x 2 ) 



x + x 3 - x 2 <9$ . 
(x , Xi, x 4 ) - -; - — (x , x 3 , x 2 ) 



(xi - x 3 )(x 2 - x 4 ) <9x ' ' (xi - x 3 )(x 2 - x 4 ) 9x 

xq + x 3 — x 4 (9$ 



(xi - x 3 )(x 2 - x 4 ) <9x 



(x ,x 3 ,x 4 )} , 



(B7) 



with A = - 7j5 + ln(4vrxj, A ± = 3 - 2 7e + 21n^. Here, x m = /i^ p /m 2 , and /i NP is the 
scale to integrate heavy particles out. 

Defining A 2 = Xq + x 2 + x 2 — 2x xx — 2x x 2 — 2xiX 2 , the two-loop vacuum function 
$(x ,Xi,x 2 ) is written as 
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A 2 > 0, y/xi + < v^o : 

&(x , Xi, x 2 ) = (x + xi - x 2 ) lnx hixi + (x — x\ + x 2 ) lnx lnx 2 
+ (— xq + sci + x 2 ) lnxi In x 2 + A|2 In 1 



2x 

/x -xi+x 2 -A\ , xi x 2 or /x + Xi-x 2 -A 

x In ( J - In — In 2L l2 ( 

v 2x ' x x v 2x 



where L i2 (x) is the spence function; 
A 2 > 0, -y/xo + t/x^ < y/x{: 

$(x ,Xi,x 2 ) = Eq.(jB8l)(xo «-> xi) ; (B9) 

A 2 > 0, v /x^ + y/x[ < ^fx~ 2 : 

$(x ,xux 2 ) =Eq.dBHI)(xo^x 2 ) ; (BIO) 

A 2 < 0: 

$(x , Xi,x 2 ) = (x + x\ — x 2 ) lnx lnxi + (x - x± + x 2 ) lnx hix 2 

+ (-x + x 1 + x 2 ) In a* lnx 2 - 2^{C7 2 ( 2 arccos( ~ X ° + Xl + * 2 



2y/XiX 2 

" C/2 ( 2arCC ° S( 2^T 2 ))+^( 2 ™( ))} • ( B11 ) 



where CU(x) denotes the Clausen function. 
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